Spin relaxation in semiconductor quantum dots 
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The spin relaxation time due to the electron-acoustic phonon scattering in GaAs quantum dots 
is studied after the exact diagonalization of the electron Hamiltonian with the spin-orbit coupling. 
Different effects such as the magnetic field, the quantum dot size and the temperature on the spin 
relaxation time are investigated in detail. Moreover, we show that the perturbation method widely 
used in the literature is inadequate in accounting for the electron structure and therefore the spin 
relaxation time. 

PACS numbers: 71.70.Ej, 73.21.La,72.25.Rb 



I. INTRODUCTION 



Spin-related phenomena in semiconductors have at- 
tracted much attention recently as they are the key in- 
gredient in the field of spintronics. 1 Among these, the 
spin-orbit coupling mechanisms in semiconductor quan- 
tum dots (QDs) provide a basis for device applications 
such as qubits in quantum computers and have there- 
fore caused much interest. 2 ' 3 ' 4,5,6 ' 7 Voskoboynikov et at 
studied the electron structures of QDs by exactly diago- 
nalizing the Hamiltonian with spin-orbit coupling. 8 Gov- 
ernale studied the electron structure of few-electron in- 
teracting QDs with Rashba spin-orbit coupling by spin 
density functional theory. 9 Valin-Rodriguez et al. in- 
vestigated spin procession in QDs with the spin-orbit 
coupling. 2 Besides the effect of the spin-mixing in the 
electron structure, the spin-orbit coupling also induces 
the spin relaxation via further coupling with phonons, 
which alone conserve the spin and therefore are unable 
to cause any spin relaxation. Many works calculated the 
spin relaxation time (SRT) due to the spin-orbit coupling 
induced spin-flip electron-phonon scattering at zero or 
very low temperatures. 4,6,7,10 Unlike the electron struc- 
ture calculation, 2,8,9 to our knowledge all works on the 
SRT are based on perturbation theory where the spin- 
orbit coupling is treated as a perturbation in the Hilbcrt 
space spanned by Hq which does not include the spin- 
orbit coupling. Moreover only the lowest few energy lev- 
els of Hq are included in the theory. 4,6,7,10 Whether the 
perturbation based on the lowest few levels of Hq is ad- 
equate remains unchecked. 

In the present paper, we investigate the SRT of GaAs 
QDs confined in the quantum well by parabolic poten- 
tials by exactly diagonalizing the total Hamiltonian. We 
calculate the SRT due to the scattering with the acous- 
tic phonons by the Fermi golden rule after getting the 
energy spectra and the wavefunctions from the exact di- 
agonalization. We find that the perturbation approach 
is inadequate in calculating the SRT and therefore it is 
necessary to reinvestigate the accurate SRT via the ex- 
act diagonalization approach. We organize the paper as 
follows: In Sec. II we set up our model and the Hamilto- 



nian. Then in Sec. Ill we present our numerical results: 
We first compare the results obtained from our exact di- 
agonalization method with those from the perturbation 
approach and show that the perturbation method is in- 
adequate in accounting for the SRT in Sec. III(A). In 
Sec. III(B) we discuss the SRT of a small QD at T = 4 K 
where only the lowest two energy levels of the total elec- 
tron Hamiltonian after the exact diagonalization account 
for the SRT. Nevertheless at least 12 energy levels of Hq 
are needed to get these lowest two energy levels. We then 
turn to effects of the magnetic field, the temperature and 
the quantum well width to the SRT in Sec. 111(C)- (E). 
We give our conclusions in Sec. IV. 



II. MODEL AND HAMILTONIAN 

We set up a simplified model to study the spin relax- 
ation in the QD's which are defined by parabolic poten- 



tials V c (r) 



\v in a quantum well of width a. A 



magnetic field B is applied along the growth (z) direction 
of the quantum well. The total Hamiltonian is given by 



H = H P + H 



ph 



H P 



(1) 



with the electron Hamiltonian H e — Ho + H so . Here Hq 
is electron Hamiltonian without the spin-orbit coupling: 



2m* 



K(r) + H B 



(2) 



in which P = -i%V + (e/c)A with A = ~(-y,x,0) 
stands for the electron momentum operator, m* is the 
electron effective mass. Hb = \q^bP^z is the Zeeman 
energy with ^representing the Pauli matrices. H so = 
•yh ■ cr is the spin-orbit coupling which is the key to the 
spin flip and spin relaxation, h = [P x {Py — P*), Py(P z ~ 
P 2 ), P z {Px ~ Py)\ i s t ne Dresshauls effective magnetic 
field in the bulk material. 11 In quantum well with small 
width, H so can be simplified as 



7c(- 



PyOy) 



(3) 



2 



with 7 C = 7(7r/a) 2 . H p h in Eq. (1) is the Hamiltonian for 

phonons and is given by H p h — X) q A ^• l - Ll qA a qA a qA with 
ujq\ standing for the phonon energy spectrum of branch 
A and momentum q. The electron-phonon scattering is 
given by 

H ?p = E M q A ( a qA + a q-0 cx P(*q ' r ) ( 4 ) 

qA 

with M q A being the scattering matrix element. 

We diagonalize the electron Hamiltonian H e in the 
Hilbert space \n, I, a) constructed by Hq = j-tt + V c (r) + 
Hb: \*t) = E n i„C e nl(T \n,l,a). Here H \n,l,a) = 
E ni i t<r \n,l,a) with 

(v\n,l,a) = N ntl (ar)We-^LW((ar) 2 )e iW X(T , (5) 
and 

E nM = m{2n + \l\ + 1) - hu B l + (tE b . (6) 
I 



In these equations n = 0, 1, 2, • • • and / = , ±1, ±2, • • • 

are quantum numbers, fl = yu>o + and ujg = 

i 

eB/(2m*). N n>l = with a = ^HTJh. 

Eb = \g^B^ is the Zeeman splitting energy, a = ±1 
refers to the spin polarization along the z-axis. \ a rep- 
resents the eigenfunction of a z . l)n is the generalized 
Laguerre polynomial. By solving 



H e \*t) = e e \V e ) , (7) 



one can determine the eigenenergy e£ and the eigenfunc- 
tion of the total electron system H e . It is noted that 
due to the presence of the spin-orbit coupling H SOl a is 
no longer a good quantum number. Mixing occurs for 
opposite spins: 



(n,l,a\H so \n',l',a') = i2n lc aS l , +(T , l 6 a ,^[a(co B /n)A^ n , ihl , - <rl' A { ^ n , I V + A^ n , I V ] . (8) 
It is this mixing that makes the originally spin-conserving e lectron-phonon scattering Eq. (4) cause spin relaxation. 



to A { V in Eq. (8) are given in detail in Appendix A. 
The eigenfunction |\f^) obtained from Eq. (7) contains 
spin mixing for each state I. We assign an eigen state 
I to be spin-up if u z = (^e\a z \^e) > or spin-down if 
a z < 0. An electron at initial electron state i with energy 
£i and a spin polarization can be scattered by the phonon 
into another state / with energy e / and the opposite spin 
polarization. The rate of such scattering can be described 
by the Fermi golden rule: 

2-7T 

qA 

+ (nqX + l)S(e f -e l +uj ciX )} , (9) 

with n q A representing the Bose distribution of phonon 
with mode A and momentum q at the temperature T. 
Its expression after the integration is given in Appendix 
B. The SRT t can therefore be determined by 

i f 

in which /j = Cexp[— e i /(fc^T)] denotes the Maxwell 
distribution of the i-th level with C being a constant. 

III. NUMERICAL RESULTS 

We perform a numerical investigation of the SRT [Eq. 
(10)] in GaAs quantum dots at low temperatures by di- 



agonalizing the Hamiltonian H e for each given dot size: 
the quantum well width a and the effective diameter 

d = m*L a 1 as we ^ as eacn gi yen applied magnetic field 
B. To do so, we gradually increase the number of basis 
function arranged in the order of energy in the Hilbert 
space \n,l,a) to ensure the 0.1 % precision of the con- 
verged energy q. As an example for a QD with B = 1 
T and a = 5 nm, when d = 20 nm, in order to converge 
the lowest 2 (100) levels, one has to use 12 (120) basis 
functions; nevertheless, when d = 60 nm, in order to con- 
verge the lowest 2 (100) levels, one has to use 20 (200) 
basis functions. 



D 


5.3 x 10 3 kg/m 3 


K 


12.9 


Vst 


2.48 x 10 3 m/s 


g 


-0.44 


V s l 


5.29 x 10 3 m/s 




7.0 eV 


ei4 


1.41 x 10 3 V/m 




0.067m 



TABLE I: Parameters used in the calculation 

The electron-phonon scattering is composed of the fol- 
lowing contributions: (i) The electron-acoustic phonon 
scattering due to the deformation potential with M^ s[ = 

2 py - ; (ii) The electron-acoustic phonon scattering due 
to the piezoelectric field for the longitudinal phonon 

mode with Ml pl = 32 K f^f 4 < 3 *y >' and for the 

two transverse phonon modes with X^j=i i ^-qpt, = 



3 



fS#[<7^ + till + fal - Here S stands 

for the acoustic deformation potential, D is the GaAs 
volume density, eu represents the piezoelectric constant 
and k denotes the static dielectric constant. The acoustic 
phonon spectra ujqx are given by uiqi = v s iq for the lon- 
gitudinal mode and w qpt = v st q for the transverse mode 
with v s i and v a t representing the corresponding sound ve- 
locities. It is noted that notwithstanding the fact that we 
include all these acoustic phonons throughout our com- 
putation, for all the cases we have studied in this paper, 
the main contribution comes from the electron-phonon 
scattering due to the piezoelectric field for the transverse 
mode with the later being at least one order of magni- 
tude larger than the other phonon modes. Moreover, the 
contribution from LO-phonon is negligible in the temper- 
ature regime we are studying. The parameters used in 
our calculation are listed in Table I. 13 



A. Comparison with previous works at T — 4 K 

We first compare our approach with the perturba- 
tion calculations widely used in the literature 6 ' 7 ' 10 at low 
temperature to double check the validity of our method 
as well as that of the perturbation method where H so 
is treated as the perturbation. Following the previous 
works, 6,7 ' 10 we calculate the SRT between the lowest two 
Zeeman splitting levels at d = 20 nm and T = 4 K. Un- 
less specified, the width of the quantum well a is fixed to 
be 5 nm throughout the paper. To the first order of H so , 
the energy difference of the lowest two states with the 
opposite spins is AE — 2Eb (the first order correction is 
zero) and the wave functions of these two states are 



*T = (r|0,0,T) , 

* x = (r|0,0a)-^<r|0,l,T) , 

^2«(l-eB/(2fia 2 )) 



(11) 



in which A = ihrf c {-) ^ ' ^ » . The SRT r is 

"-VO' iio,l,T — ^0,0,1 

therefore given by 



- = c\A\ 2 n q q 6 / d6>sin b <9(sin 4 (9 + 8cos 4 (9) 



x exp(— -q 2 sin 2 9)I 2 (qcost 



(12) 



with q = AE/(hv st a), c = 9irae 2 e 2 4 /(HDv 2 t K 2 ) and 
I(q z ) = 8Tr 2 sin(aq z /2)/{aq z [^ 2 - (a^) 2 )]}- This is ex- 
actly the same calculation used in the literature. 6 ' 7 ' 10 

In Fig. 1 we compare the SRT calculated from Eq. (12) 
(curve with ■) with our exact diagonalization method 
(curve with O) described in the previous section, but 
with only the lowest four levels of the Hilbert space, 
ie., |0,0, a) and |0, 1,<t) (<r =| or J.), taken, correspond- 
ing to the same levels used in the perturbation method. 
It is seen from the figure that there is at least one or- 
der of magnitude difference between the two curves: r 
obtained from the perturbation method is much larger 
than the one from the exact diagonalization method. 



Moreover, the trends of the magnetic field dependence 
are also different. We point out that these differences 
arise from the fact that only the first order of the per- 
turbation is applied. It can be fixed if one further in- 
cludes the second order correction in AE, ie. AE = 
2E B + \A\ 2 (E ,i^ — £0,0,4.)- 1* 1S noted that the second 
order correction to the energy difference is much larger 
than 2Eb, the unperturbed energy difference. The SRT 
calculated from the perturbation method, modified with 
the energy correction to the second order, is plotted as a 
function of the magnetic field B in the same figure (curve 
with x). One notices that it exactly hits on the curve 
from our diagonalization method (curve with Q)- 
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FIG. 1: SRT versus the magnetic field. Curve with ■: per- 
turbation result without energy corrections; Curve with : 
diagonalization result but with only the lowest four levels used 
as basis functions; Curve with x : Perturbation result with the 
second order energy correction;Curve with ♦: exact diagonal- 
ization result with the energy sufficiently converged; Curve 
with □: diagonalization result with the lowest six levels used 
as basis functions; Curve with H 4 : perturbation result with 
the lowest six levels of Ho as basis functions and with the 
second order energy correction. 

We notice that the diagonalization above includes only 
the lowest four levels. If it is adequate in converging the 
lowest two levels remains unchecked. As mentioned in 
the beginning of this section, that for the size of the dot 
we are studying here, in order to converge the lowest two 
energy levels, one has to use 12 basis functions, r cal- 
culated from the exact diagonalization method is plotted 
in also Fig. 1 against the magnetic field (curve with ♦). 
Strikingly, it is orders of magnitude larger than that from 
the perturbation. 

In order to understand this huge difference, now we in- 
clude six lowest energy levels of H , i.e., |0, 0, |), |0, 0, J.), 
10,1,4), |0, 1,|), |0, — 1,T) and |0, -1,4), as basis func- 
tions in the perturbation method. The wavefunctions of 
the lowest two states of H e are therefore given by 



*T = (r|0,0,T) 
*4 = MO, 0,|) 



■B<r|0,-U) , 
-4<r|0,l,T) , 



(13) 



4 



in which 



B = ihcrf c (-) 



71\ 2 1 



eB/(2ha 2 



(14) 



' « ^o,-i,T _ ^0,04 
The energy difference between vpf and ^ now becomes 

A.E = 2Eb + \A\ 2 (Eoa^— Eofi,i) — |^| 2 (-Eo : -i,t ~ ^o,o,t) 
The corresponding SRT r is hence given by 



(15) 



1 

T 



c\A-B\ 2 n, 



Jo 



d6>sin 5 6>(sin 4 (9 + 8cos 4 (9) 



x cxp(--q 2 sin 2 6)I 2 (q cos 6) , 



(16) 



The numerical results of Eq. (16) are plotted by the curve 
with in Fig. 1 . It is seen from the figure that the inclu- 
sion of the additional basis functions in the perturbation 
method also greatly enhances the SRT by orders of mag- 
nitude. The results obtained from the exact diagonaliza- 
tion method with the same lowest six levels as basis are 
given by the curve with □. The two curves are almost 
the same, and much closer to the final converged results 
(curve marked with ♦). 
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FIG. 2: SRT versus the applied magnetic field at d=20 nm 
and T=4 K. The inset is the corresponding energy splitting 
AE between the lowest two levels 



lowest two levels (for a QD with d — 20 nm at 4 K, the 
SRT is determined by the lowest two levels of the total 
electron Hamiltonian H e ), one has to use a lot of basis 
functions to converge the energy and the resulting SRT 
between these two levels is therefore strongly readjusted. 
This is because the spin-orbit coupling is very strong in 
mixing different energy levels of Hq. 

In the following subsections, we therefore reinvestigate 
the properties of the SRT based on the exact diagonal- 
ization calculation. 



B. SRT of a d = 20 nm QD at T = 4 K 

As pointed out in the previous subsection that for a 
QD with d = 20 nm, at T = 4 K the SRT is determined 
by the spin-flip transition between the lowest two energy 
levels after the exact diagonalization, although at least 
12 energy levels in the Hilbert space of Ho are essential 
in getting these two levels. In this subsection we focus on 
the effects of the external fields on the SRT determined 
by these two levels. 

It is seen from Eqs. (12) and (16) that in a given basis 
the spin relaxation rate 1/r is determined by two com- 
peting trends as a function of the energy splitting AE: 
(i) q n q , which increases with AE in the present case, 
and (ii) exp(— ^q 2 sin 2 8)I 2 (qcos0), which decreases with 
AE. Therefore, the SRT can be uniquely determined by 
the energy AE. For small AE, it is easy to see that the 
trend (i) dominates the when AE < 7.0hv st /d, which is 
0.57 meV at d = 20 nm. That is, the SRT decreases with 
AE when AE < 0.57 meV. 

In Fig. 2 the SRT is plotted against the applied mag- 
netic field B. It is seen from the figure that r decreases 
with the applied magnetic field. This is understood from 
the fact that the energy splitting AE increases with the 
applied magnetic field as shown in the inset. Moreover, 
even for the largest energy splitting 0.04 meV at B = 1 
T, it is one order of magnitude smaller than 0.57 meV, 
energy splitting required to have the opposite t-B de- 
pendence. 



C. Magnetic field dependence of the SRT 



It is clearly seen from the above calculation that the 
perturbation approach widely used in the literature is 
inadequate in describing the SRT even with the second 
order energy corrections included. In principal in order 
to use the perturbation method to calculate the SRT, 
one has to include sufficient number of the states in the 
basis in stead of only the lowest four levels widely used 
in the literature. This is of course inapplicable especially 
for larger QD's or higher temperature where one has to 
include a lot of basis functions (for a QD of d = 60 nm, 
one has to use 100 levels as basis functions) and the SRT 
is determined by many levels (in stead of only the lowest 
two) of the total electron Hamiltonian H e . Even for the 



We investigate the magnetic field dependence of the 
SRT for different diameters of the QD's at two different 
temperatures as shown in Fig. 3. Unlike the previous 
subsection where only the lowest two energy levels are 
important, here for most cases one has to include many 
levels of the total electron Hamiltonian. 

It is seen that the SRT decreases rapidly with the mag- 
netic field at each dot size and temperature. This feature 
is quite opposite to the bulk, 14 the two- 15 and the one- 
dimensional 16 cases where the SRT always increases with 
the magnetic field. This is because in the dot case there 
are only discrete energy levels and the magnetic field 
helps to increase the spin-flip scattering as discussed in 
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FIG. 3: The SFT vs. the magnetic field with different sizes of 
the quantum dots:(B) d= 20 nm,(») d= 30 nm, (a) d= 40 nm, 
and (♦) d= 60 nm for T = 4 K (a) and 10 K (b). 
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FIG. 4: The SRT vs. the temperature under different mag- 
netic fields at d = 40 nm. Curve with ■: B = 0.1 T; Curve 
with •: B = 0.5 T; and curve with A: B = 1 T. 



temperatures, the occupation to the high energy levels 
becomes important and it is inadequate to consider only 
the lowest several levels. For lower magnetic fields, the 
space between different energy levels is smaller. There- 
fore, more levels are included in the energy regime deter- 
mined by fi in Eq. (10) which leads to a faster response 
to the temperature. This feature is more pronounced 
in the low temperature regime. For high temperatures, 
as there are already many levels included in the energy 
space, adding a few more levels does not change the SRT 
significantly. Consequently the rates of the decrease of 
the SRT with the temperature become similar for differ- 
ent magnetic fields when T > 16 K. 



the previous subsection. Moreover, one notices that the 
SRT drops dramatically with the dot size. For a dot with 
d = 60 nm, the SRT is more than 6 orders of magnitude 
faster than the one with d = 20 nm. This is understood 
that for larger dots, more energy levels are engaged in the 
spin- flip scattering and hence sharply reduce the SRT. 

D. Temperature dependence of the SRT 

We plot the SRT as a function of the temperature in 
Fig. 4 for a QD with d — 40 nm under three different 
magnetic fields. From the figure one finds that the SRT 
gets smaller with the increase of the temperature. More- 
over, the smaller the magnetic field is, the faster the SRT 
drops with the temperature. 

These features can be understood as follows: With 
the increase of the temperature, the phonon number n q A 
gets larger. This enhances the electron-phonon scatter- 
ing and leads to the larger transition probability. More- 
over, unlike the previous work 10 where the difference be- 
tween zero temperature and finite temperatures is just 
the phonon Bose distribution, we stress that for high 



E. Well width dependence of the SRT 
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FIG. 5: The SRT vs. the width of the quantum well in differ- 
ent temperatures at B= 1 T. Curve with ■: T = 4 K; Curve 
with •: T = 10 K; and Curve with A: T = 20 K. 
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As the QDs are confined in the quantum well, it is 
necessary to study the quantum well width dependence 
of the SRT as shown in Fig. 5 where r is plotted as a 
function of the well width a for different temperatures 
at B=l T. It is noted that the SRT increases with the 
well width a. This is due to the fact that the spin-orbit 
coupling H so [Eq. (3)] is proportional to 1/a 2 . Smaller 
well width corresponds to larger spin-orbit coupling and 
therefore smaller SRT. We point out here that the well 
width in the present calculation is much smaller than the 
dot size d and only the lowest subband contributes to the 
SRT. For larger well width, more subbands are involved 
and hence there adds an opposite tendency for a shorter 
SRT with the increase of the well width. 



IV. CONCLUSIONS 

In conclusion, we have investigated the SRT in GaAs 
QDs by the exact diagonalization method with applied 
magnetic fields. After comparing the exact diagonaliza- 
tion method with the perturbation approach widely used 
in the literature, we find that the later is inadequate in ac- 
counting for the electron structure and the SRT in QDs. 
This is because that the energy splitting caused by the 
spin-orbit coupling is several times large than the Zee- 
man splitting used in the perturbation approach. More- 
over, a lot more energy levels of Hq are coupled by the 
spin-orbit coupling and therefore contribute to the lowest 
energy levels of the total QD Hamiltonian. We therefore 
reinvestigated the SRT from the exact diagonalization 
method to explore its dependence on the magnetic filed, 
the temperature and the size of the QD. We find the SRT 
decreases with the magnetic field, which is quite opposite 
to the bulk, the two- and one- dimensional cases. It also 
decreases with the diameter of the QD, but increases with 
the width of the quantum well on which the QD grows. 
For high temperature, the SRT becomes much faster due 
to the stronger electron-phonon scattering and the wider 
range of energy space the electron occupies. All our inves- 
tigation suggests the importance of the exact calculation 



of the energy structure. 
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APPENDIX A: THE EXPRESSIONS OF A w , A (2) 

AND A (3) 

and A^ in Eq. (8) are given by 

/>oo 

<Uf = oc \ r 2 R n Ar)Rn'M(r)dr , (Al) 
Jo 
1 f°° 

Ox.f = - / Rn,i{r)Rn>Ar)dr , (A2) 
a Jo 

<!,™m' = -/ o rR n , l (r)-R n/ , l/ (r)dr , (A3) 

where R n j = ^^(or)!'! cxp (-i^l)L^(a 2 r 2 ) is 

the spatial part of the wave function Eq. (5). From the 
integration over the angular part, we get the relation |Z — 
l'\ = 1. Substituting this relation into Eqs. (Al)-(A3), 
after carrying out the integration we have 

A n},n>,l> = ^-Wn+\l\ + lS n , n > - >M,',n-l) , (A4) 



^0, otherwisi 



It is noted that due to the symmetry between {n, 1} and {n', I'}, in above two equations we only give the results 
with \l'\ = \l\ + 1. Finally 



,0) 



-l/'l/t (2) 4- A (1> -i 



(1) 



( Vn+m 
in 



n'\(n+\l\)\ 
n!(n' + |Z'|)! - 



lo, 



if |Z'| = \l\ - 1 and n' = n 

if |Z'| = \l\ + 1 and n' > n 

if |Z'| = |/| + 1 and n' = n- 1 
otherwise 



(A6) 
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APPENDIX B: THE EXPRESSION OF 



with q = (Qcos6>, Qsind, q z ) and g = |q| = 1 ' ■ Here AT, = n 9 if ^ > £7, or n q + 1 if ^ < £j. Gjj in Eq. (Bl) 
is 

G iJ {Q 2 /{Aa 2 ),q z )= ]T ^^(C^)* <«2, k| exp [t( fea: + Zx> exp , (B2) 



k 2a 2/ 

ni,ii,n2,(2,CT 



in which 



Q2 

exp(^2)(n 2 ,/ 2 |exp[i(g x a; + Q'j / y)]|ni,Zi) = (B3) 

/ i i n Ui— (2I ni ™ 2 /o2 

V(ni + |/i|)!(n2 + N)! ^ l 1 " a ' |Ia| " 

with sgn(a;) denoting the sign function, C™ ; = - (^Z^) and n = i + j + 

I ' 
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